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ABSTRACT 

Aims. The aim of the project is to search for lithium in absorption at 6707.8A to constrain the nature and the mass of the brightest 
low-metallicity L-type dwarfs (refered to as L subdwarfs) identified in large-scale surveys. 

Methods. We obtained low- to intermediate-resolution (R ~ 2500-9000) optical (~560-770 nm) spectra of two mid-L subdwarfs, 
SDSSJ125637.13-022452.4 (SDSS1256; sdL3.5) and 2MASS J162620.14-H392519.5 (2MASS1626; sdL4) with spectrographs on 
the European Southern Observatory Very Large Telescope and the Gran Telescopic de Canarias. 

Results. We report the presence of a feature at the nominal position of the lithium absorption doublet at 6707.8A in the spectrum 
of SDSS1256, with an equivalent width of 66±27A at 2.4cr, which we identify as arising from a CaH molecular transition based on 
atmosphere models. We do not see any feature at the position of the lithium feature in the spectrum of 2MASS1626. The existence 
of overlapping molecular absorption sets a confusion detection limit of [Li/H] = -3 for equivalently-typed L subdwarfs. We provided 
improved radial velocity measurements of -126±10 km s“' and -239±12 km s“' for SDSS1256 and 2MASS1626, respectively, as 
well as revised Galactic orbits. We implemented adjusting factors for the CaH molecule in combination with the NextGen atmosphere 
models to fit the optical spectrum of SDSS1256 in the 6200-7300A range. We also estimate the expected Li abundance from inter¬ 
stellar accretion ([Li/H] = -5), place limits on circumstellar accretion (10® g/yr), and discuss the prospects of Li searches in cooler L 
and T subdwarfs. 

Key words. Stars: brown dwarfs — Stars: subdwarfs — techniques: spectroscopic 


1. Introduction 


Stars spend most of their lifetime on the main-sequence in hy¬ 
drostatic equilibrium, burning hydrogen. To the contrary, brown 
dwarfs never reach c ore temperatur e s and pressures high enoug h 
to fuse hydrogen (iKuman Il963t iHavashi & Nakanol Il96a) . 
Lithium (Li) is destroyed in stellar interiors via collisions with 
protons. At solar metallicities brown dwarfs with masses below 
0.065 M(7) do not reach the temperatures needed for Li I bu rning 

dMagazzu et alJl9^lRebolo et al.lfT992l : lBasri et alJl996l) . The 

precise mass below which Li I is preserved depends on metallic- 
ity and, according to models , it cou l d be higher at lower metal¬ 


licities (lBurrows & Liebig I1993I: Qiabrie£&^Muff^Jl^97t 
Bildsten et al. 1997t iBaraffe et alJ 1998I : Ichabrier & Baraffel 
^OOUBurke et alJl2004r 


Low-resolution optical spectroscopy of solar-metallicity 
nearby L dwarfs suggests early-L dwarfs represent a mixture 
of low-mass stars and brown dwarfs whereas most objects with 
spectral types later than ~L5 are substellar. F rom a large sam- 
ple of nearby L dwarfs with optical spectra, iKirknatrick et alJ 


* Based on observations collected at the European Southern 
Observatory, Chile, under programmes 089.C-0883 and 09LC-0594A 
** Based on observations made with the Gran Telescopio Canarias 
(GTC), installed in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias, in the island of 
La Palma (programmes GTC64_10B and GTC38_1 lA) 


(l2008h showed that 10-20% of L0-L2 dwarfs exhibit Li in 
absorption, the percentage increasing to 40-80% for mid- to 
late-L dwarfs. This trend is corr o borated by o t her works, e.g. 


2000I). ^vlenkoetd] 


iTinney (1 19981). IKirknatrick et al.l( 199 91 1 
(l2007a ). ICruz et alT ( 20091) . and IZapatero Osorio et al.l ( 20141) . 
In the T dwarf regime, only Luhmanl6B is known to har¬ 
bour Li in a bsorpt i on with an equiva lent width of ~8A 


(iFahertv et al.l 1201 4 lLodieuet^lJ_^0]^ whereas Li „ _ 

detec ted in G1229 B (ISchultz et al.l 19981 : lOnnenheimer et al 
1998 ) and in the two component of Epsilon Indi B dKing et al 


20101) . Lithium remai ned undetected in the optical spectrum of 


another 17 T dwar fs (iBurgasser et al.ll^03bl: iKirkpatrick et al.l 
iTOOilLeggett et al.ll2012l). 

Ultracool subdwarfs (or Population II ultracool dwarfs) 
are defined as metal-poor dwarfs with spectral types later 
than M7. They appear bluer than the main-sequence of solar- 
metallici ty stars due to the dearth of metals in their atmo¬ 
spheres (iMould & Hylandll 19761: IBaraffe et al.lll997l) . They usu¬ 
ally exhibit halo kinematics, i ncluding hig h proper motions and 
large heliocentric yelocities (lGizisl[l997l) . They are important 
tracers of the chemical enrichment history of the Galaxy be¬ 
cause they belong to the first generations of stars. The cur¬ 
rent compendium of ultracool subdwarfs contain s approximately 


50 obj e cts, in part i cular 1 1 L-type subdwarfs (IBurgasser et al 
12003 ai iBurgasseil l2004 ICush^ing et al.l 120091 : Siyarani et al 
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Table 1. Properties of the tw o mid-L subdwarfs present ed in this paper. Listed are the names, coordinates, spectral types ([Burgasseil 
l2007t iBurgasser et aL 2009ll. Sloan r-band magnitude dAdelman-McCarthv et al.ir2007t lAdelman-McCarthv & et al]l2009h . Droper 
motions, distances dSchilbach et alj|2009h . final upper limits on the Li I pseudo-equivalent widths, radial velocities, and revised 
space motion (this paper). 


R.A. dec 

SpT 

SDSSr 

fiaCOSd 


dist 

EWu 

RV 

(U,V,W) 

hh:mm:ss.ss 


mag 

mas/yr 

mas/yr 

pc 

mA 

km s“' 

(km/s,km/s,km/s) 

12:56:37.13 -02:24:52.4 

sdL3.5 

21.82 

-741.1±1.4 

-1002.0±1.4 

53.3±6.8 

<66±27 

-126±10 

(-62±6,-227±37,-224± 18) 

16:26:20.14 -1-39:25:19.5 

sdL4.0 

20.65 

-1374.1±1.0 

-h238.0±0.9 

33.5±1.3 

<90 

-239±12 

(-166±6,-261±9,-2±11) 


2009 

; Lodieu et al.l 

2010]; 

Kirknatrick et al.l 120101: iLodieu et al.l 

2012 

;iKirknatrick et al.||2014). 


Because M subdwarfs are old components of our Galaxy (a 
few Gyr) with effective temperatures in the 3000^000 K range 
(IRaipurohit et al.l 120141) we do not expect to detect lithium in 
absorption in their optical spectra because all the Li has been 
burnt. However, L subdwarfs are very low mass objects at low 
metallicities and are likely to have masses close to the minimum 
needed for Li I burning. Lithium is expected to be present in 
the formation of subdwarfs at high abundances given its origin 
in Big Bang nucleosynthesis (ISpite & Spitelll98^ iRebolo et al.l 


Il988t [Bonifacio et alJl200% . Hence, L subdwarfs might have 
preserved a significant amount of primordial Li I in their at¬ 
mospheres if they are substellar. To enable the detection of the 
Li I resonance doublet at 6707.8A, it is advisable to focus on 
the coolest subdwarfs known, with effective temperatures be¬ 
low ~2000-2500K. However, the faint continua of L subdwarfs 
around the Li I feature at 6707.8A necessitates observations of 
sufficient resolution and sensitivity to firmly establish its pres¬ 
ence or absence. For depletion factors in the range 2-100, theo¬ 
retical models predict equivalent widths (EW) between 0.3A and 
0.8A (iPavlenko et al.|[T99^ . Prior, low-resolution optical obser- 
va tions did not have suffi c ient signal-to-noise t o do this (Table 
2l BumassereLak 2003al: ICushing et al.ll2009t IBurgasser et al.l 
2009HKirkDatrick et al.ll2010lK 


In this paper we present medium-resolution optical spectra of 
two mid-L subdwa rfs, SPSS J125637.13-022452.4 (hereafter 
SDSS1256; sdL3.5: LS i varan let al.ll2009t IBurgasser et al.ll2009h 
and 2 MASSJ162620 1 4-H39 2519.5 (hereafter 2MASS1626; 
sdL4; iBurgasser et al.l l2007h . obtained with the European 
Southern Observatory (ESO) Very Large Telescope (VLT) in 
Paranal (Chile) and the Gran Telescopio de Canarias (GTC) in 
the Roque de Los Muchachos Observatory on La Palma (Canary 
Islands). We selected these two subdwarfs because they repre¬ 
sent the brightest of their class. Their properties are summarised 
in Table [1] In Section |2] we describe the spectroscopic observa¬ 
tions conducted in service mode by ESO and GTC. In Sections 
[3 and |6] we present refined radial velocities yielding improved 
3D space motions as well as luminosities and effective tempera¬ 
tures for both objects, respectively. In Section!?] we discuss our 
analysis of the lithium absorption in the optical spectra of both 
mid-L subdwarfs. In Section|3 we compare our upper limits on 
the pseudo-equivalent widths of the Li I feature with theoreti¬ 
cal predictions to place a limit on the mass of the subdwarfs. In 
SectionQwe place our results in context and discuss the oppor¬ 
tunity offered to detect lithium in halo brown dwarfs. 


2. Spectroscopic observations 

2.1. GTC/OSIRIS optical spectroscopy 


We obtained optical spectroscopy of 2MASS1626 with the 
OSIRIS (Optical System for Imag ing and low-inte rmediate 
Resolution Integrated Spectroscopy; ICena et ^120001) spectro¬ 
graph mounted on the 10.4-m GTC telescope in La Palma. 
Observations were conducted in February and April 2011 over 
two semesters under programmes GTC64_10B and GTC38_1 lA 
(PI Espinoza Contreras; TableO. 

OSIRIS is equipped with two 2048x4096 Marconi CCD42- 
82 detectors offering a field-of-view approximately 7x7 arcmin 
with an unbinned pixel scale of 0.125 arcsec. We employed the 
R2500R grating with a slit width of one arcsec giving a reso¬ 
lution of ~2400 under dark or grey conditions, variable seeing 
between 0.8 and 1.6 arcsec, and clear skies. We obtained three 
exposures shifted along the slit with on-source integrations of 
1402 s on 5 February 2011, another three with the same on- 
source integrations on 12 February 2011, and two additional 
exposures of 1636 s on 13 April 2011. Bias frames, dome flat 
fields, and Ne, Xe, and HgAr arc lamps were observed by the 
observatory staff during the afternoon preceding or following the 
observations. 

We carried out the data reduction of the o ptical spectra un¬ 
der the IRAF environment (lTodvlll9^ Il993ln . We subtracted 
the raw median-combined bias frame from the raw spectrum and 
divided by a normalised dome flat field. We extracted each spec¬ 
trum individually without combining them to avoid shifts due to 
the large space velocity of the target and the position of the Earth 
at the time of the observations. We chose an optimal sky back¬ 
ground and aperture to extract the one-dimensional (ID) spec¬ 
trum. We calibrated our spectrum in wavelength with arc lamps 
with an accuracy better than 0.05A rms. Einally, we calibrated 
each ID spectrum in flux with the spectro-photometric standard 
stars G191-B2B (DA . 8; Ivan LeeuwenI 1200^ iHpg etaP 120001: 
Gianninas_e^^ 1201 Ih a nd Ross 640 (DZ5: [Harrington & Dahra 


19801 : ISion et aP 120091 : iLenine & Sharal 1200.51) . The individual 


spectra were shifted by the barycentric velocities correspond¬ 
ing to the times of observations, then combined with uncertainty 
weighting. The final combined optical spectrum of 2MASS1626 
is displayed on the right-hand side panel in Eigure[T] The signal- 
to-noise around the Li doublet is 90-100. 


2.2. VLT/FORS2 optical spectroscopy 

We obtained intermediate-resolution (R ~ 9000) optical spec¬ 
troscopy of SDSS1256 with the visual and near UV 


’ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation 
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Fig. 1. Full combined optical spectra of SDSS1256 (left) and 2MASS1626 (right) obtained with VLT/FORS2 and GTC/OSIRJS, 
respectively. Both spectra are in the Earth frame and have been arbitrarily normalised at 7050A. Some important bands and features 
seen in metal-poor dwarfs are marked in red. 


Table 2. Li I EWs (3cr upper limits) from low-resolution 
spectra published in the literat ure, as measured here (see 
Section B. References are: BOS ( Burgasseretalj^ 2003ah . B07 


( Burgas^r et al.l 200^ C09 (ICushing et al 
( Kirkpatrick et al .1120101) . 


20091), and KIO 


Source 

SpType 

EW(Lil) 

Reference 

A 

SDSSJ1256-0122 

sdL3.5 

<1.1 

B07 

2MASS J1626-h3925 

sdL4.0 

<0.4 

B07 

2MASSJ0616-6407 

sdL5.0 

<2.0 

C09 

2MASS J0532-H8246 

sdL7.0 

<0.6 

B03 

2MASS J0645-6646 

sdLS.O 

<3.6 

KIO 


Table 3. Log of the spectroscopic observations for SDSS1256 
and 2MASS1626 with VLT/EORS2 (1200R-H93) and 
GTC/OSIRIS (R2500) , respectively. We give the dates of 
observations and the integration times with the number of indi¬ 
vidual exposures (in case of multiple ones) for each spectrum 
obtained for each target. 

“ This spectrum is not used in our analysis due to the presence 
of a strong cosmic ray at the position of the Li I. 


Name 

Tel. Inst. 

Date 

ExpTime 



dd/mm/yy 

sec 

SDSS 1256 

ESO VLT FORS2 

27/04/12 

2824 



21/06/12 

2824 



03/04/13 

2824" 



11/04/13 

2x2824 



12/04/13 

2824 



14/04/13 

2x2824 

2MASS1626 

GTC OSIRIS 

05/02/11 

3x1682 



12/02/11 

3x1682 



13/04/11 

2x1682 


LOcal Reducer and l ow dispersion Spectrograph (LORS2; 
lAppenzeller et al.ll 19981) mounted on the ESO 8.2-m VLT at the 
Cerro Parana! observatory in Chile. We collected eight optical 
spectra over two semesters under programmes 089.C-0883A and 


091.C-0594A (PI Espinoza Contreras; Table|3]l under grey time, 
clear skies, and a seeing better than one arcsec. 

EORS2 is a multi-mode instrument mounted on the UTl 
Cassegrain focus working at optical wavelengths. We used the 
standard resolution collimator providing a pixel scale of 0.25 
arcsec and a field-of-view of 6.8 arcmin by 6.8 arcmin. We em¬ 
ployed the grism GRIS_1200R-i-93 associated to its order sepa¬ 
ration hlter GG435 with a slit width of one arcsec to achieve a 
spectral resolving power of ~9000 around the Li I doublet. We 
obtained eight exposures of 2824 s in April/June 2012 and April 
2013 (TableO, yielding a total on-source integration of ~6.3h 
and a signal-to-noise of 50-60 around 6700A. However, we do 
not use the spectrum taken on 03 April 2013 because the Li I re¬ 
gion is affected by a strong cosmic ray. Bias frames, dome flat 
helds, and arc lamps were taken by the observatory staff during 
the afternoon. 

We performed the data reduction of the LORS2 datasets in 
a similar way as the OSIRIS under IRAL. We subtracted the 
raw median-combined bias frame from each raw spectrum, di¬ 
vided each of them by a normalised dome flat field, and ex¬ 
tracted a ID spectrum calibrated in wavelength to an accuracy 
better than 0.05A rms. Linally, we calibrated each spectrum 
in f lux with the spectro-photometric stan dard stars LTT 6248 
(A: lHe>g et af] 120001: IPancino et alJ ^0121) and LTT 7379 (GO; 
van Altena et al.l 119951 ; lGontcharovl]2006t Ivan LeeuwenI 1200'^ 
Pancino et alJ |2012|) . The hnal combined optical spectrum of 


SDSS1256 is shown on the left-hand side panel in Pigure[T] 


3. Improved radial velocities and space motions 

We took advantage of the medium resolution data for SDSS1256 
and 2MASS1626 to improve their radial velocities. We hrst mea¬ 
sured the offsets of the strong Cal line at 6572.78^ averaging 
measurements from each individual spectrum. Lor SDSS 1256 
we obtained a mean offset of 1.5 (;q 3A after rejecting the value 
on 22 June 2012, which is discrepant by her. Assuming a spec¬ 
tral resolution of 0.73A equivalent to 33 km s'* for the PORS2 

^ Air wavelength from the NIST database 
http://physics.nist.gov/PhysRefData/ASD/lines_form.html 
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Velocity Shift (km/s) 


Fig. 2. Zoom on the combined spectra (black line) for SDSS1256 (left) and 2MASS1626 (right) around the Li I line. The grey 
bands indicate the formal spectral uncertainties while the horizontal error bar indicates the 3cr radial velocity uncertainties. Dotted 
lines indicate the velocity resolution of the instruments used. Overplotted at the top are marked the main features of CaH, Li, Fe, 
and TiO. We note that TiO is blanketing the entire region with cross-section of various strengths from 5x10 1x10“^®, 5x10“^® 

cm^/molecules shown in three shades of grey (light to dark). 


configuration used, we derive a mean radial velocity of -152^24 
km s“*, consistent with the values of -98. 4 km s~*, -90+4 0 
km s“', and -130+11 km s ' derived b y IWest et al.l (l2008l) . 
ISivarani et al.l (l2009t) . and iBurgasser et al.l (l2009ll . respectively. 
For 2MASS1626 we derived a mean offset of 
which translates into a radial velocity shift of - 211 . 6^33 j km s“* 
assuming a spectral resolution of L028A for the GTC/OSIRIS 
conhguration. This estimate ag rees within L5 <t with the value 
of -260+35 km s ' reported bv lBurgassed (l2004l) . 

We have also measured the radial shifts by cross-correlating 
both spectra around the 7000 -7150A TiO band to ze ro-velocity 
using SDSS templates from iBochanski et al.l (l2007l) . We used 
M7, M 8 , and M9 templates as these provided the closest 
match es to the enhanced T iO band present in these L subd¬ 
warfs (iBurgasser et al.l[2007l see Section|5]l. By correlating each 
of the individual spectra and averaging, we infer radial veloc¬ 
ities of -122+9 km s^' for SDSS1256 and -243+13 km s^' 
for 2MASS1626, where the uncertainties account for scatter be¬ 
tween observations and templates. These velocities are in agree¬ 
ment with our line measurements, but with considerably im¬ 
proved accuracy. We adopt as our hnal values the uncertainty- 
weighted mean of our line and cross-correlation measurements, 
which are -126+10 km s ' for SDSS1256 and -239+12 km s ' 
for 2MASS1626. We also verihed that these velocities aligned 
the two spectra at their rest wavelengths (Fig.|2l by checking the 
positions of a few lines (e.g. CaH at around 6698, 6705, 6717, 
and 6718A, Fel at 6710A, and TiO at 6695, 6704.5, and 6710A) 
and observed structures like to big central hump and the little 
bumps at -530 km/s and +450 km/s. 

Moreover, we measured the pseudo-equivalent width (pEW) 
of the Cal line at 6571. lOA, inferring a mean value of 3 . 54 /;Q 2 gA 
for SDSS 1256 and 3.52+® “a for 2MASS1626, in excellent 
agreement, which is expected due to the similarity in spectral 
type between both targets. 

We calculated UVW velocities for both sources based on 
our revised velocities a nd reported proper motions and distances 
(ISchilbach et al.ll2009l) . assuming a local right handed Cartesian 
coordinate system with U pointed radially outward from the 


Galactic center, V in the direction of Galactic rotation, and 
W in the direction of the Galactic North pole. We assumed 
a Local Sta ndard of Rest (LSR) v elocities of (11,12.24,7.25) 
km s ' from ISchonrich et'aH (l2010l) . Our values are reported in 
Table [H For bo t h sou r ces va lues differ from those reported in 
iBurgasser et al. (l2007l 200^ primarily due to the inclusion of 
the Schilbach et al.l (l2009l) astrometry. Both exhibit V velocities 
that are almost perfec tly counter Galac tic rotation at the Solar 
radius, 218+6 km s ' (iBovv et al.ll2012h . 

We calculated revised Galactic orbits for both sources 
using these velocitie s, following the methods outlined in 
iBurgasser et al.l(l2009l) : these are displayed in Fig. [3] SDSS 1256, 
whose galactocentric azimuthal velocity is nearly zero, exhibits 
a remarkable polar orbit that takes it completely around the 
Galactic bulge. 2MASS1626, which has a slightly retrograde az¬ 
imuthal motion and very low W velocity, remains conhned to 
the Galactic plane, but has an eccentric orbit that takes it within 
the inner kpc of the Galaxy. Neither source goes beyond a radial 
distance of 12 kpc from the Galactic center. These orbital proper¬ 
ties are consistent with membership to the inner ha lo population, 
imply ing a metallicity around [M/H] = -1.5 dex (ICarollo et al.l 
lIOOTh . 


4. The Li resonance doublet in mid-L subdwarfs 

Focusing on the Li I region, we hnd that the combined spectrum 
of SDSS 1256 shows a distinct but weak feature near the rest 
wavelength of the 6707.83A Lil line, centered within the radial 
velocity uncertainty. Using a linear fit to the adjacent peaks to 
estimate the pseudo-continuum, we measure EW = 66+27 mA 
for this feature, the uncertainty determined from the spectral er¬ 
rors and radial velocity uncertainty through Monte Carlo simu¬ 
lation. This 2.4cr detection is comparable in strength to several 
other features near the Li I line, so to assess both reliability and 
contamination we ompared this spectrum to synthetic models, as 
discussed in Section |5] 

Eor 2MASS1626, no distinct features are present, largely due 
to the much lower resolution of the data. To quantify our limits, 
we inserted a Lil absorption feature at progressively stronger 
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Fig. 3. Galactic orbits for SDSS1256 (top) and 2MASS1626 
(bottom) over ~500 Myr as viewed from above the Galactic 
North pole (left) and in cylindrical coordinates (right), fu¬ 
ture motion is indicated by solid lines, past motion by dashed 
lines, and points indicate 50 Myr steps. The Sun is located at 
(X,Y) = (-8.5,0) kpc. Representations of the Galactic bar (dark¬ 
est gray), bulge (gray), and thin disk/major spiral arms (lightest 
gray) are also shown. 



Relative Velocity (km/s) 

Fig. 4. Zoom on the spectral region containing the Li I feature 
(black), to which we added Gaussian line prohle (red) with EWs 
of 100, 200, and 400 mA, respectively. 


EWs using a Gaussian line profile modeled from the 6572.78A 
Cal line. Eigure 0] demonstrates that statistically signihcant de¬ 
viations from the original spectrum occur for EW = 90 mA, but 
a “clear” feature can only be seen for an EW > 200 mA. We are 
clearly limited by systematics associated with the low resolution 
of these data and the underlying structure of the pseudocontin¬ 
uum. The inferred limit is consistent with the best limits from 
prior low resolution observations of L subdwarfs (Table |2]i, and 
well above the strength of the feature we infer for SDSS1256. 


5. Comparison with theoreticai modeis 


To validate the possible presence of Li I in the spectrum of 
SDSS1256 and infer its atmospheric composition, we gener¬ 
ated gri ds of synthetic spectra following the procedure out¬ 


lined in iPavlenko et al 


( 2007 ah. W e starte d with the NextGen 


model atmospheres of I Allard et al.l (1200 Ih at solar metallicity. 
The prohles of the Nal and Kl resonance doublets were com- 
pu ted in the fram e work o f a quasi-static approach described 
in IPavlenko et al. (l2007bh with an upgraded approach from 
iBurrows et al. ( 2003h . Tdie line lists of VO and CaH were 
taken from Kur ucz’s websit^l with more details presented in 


Pavlenk^j20H ). The Cr H and EeH line lists w ere computed by 
Burrows et al.l (1200^ and iDulic k et al. ( 2003h . respectively. We 


upgraded the TiO line lists of iPlezI ( 1998ll with the new version 
available on his websit^B The spectroscopic data for atomic ab- 
sorption come from the Vienna Atomic Line Database (VALD; 
iKupka et al.lll999|!p i More details on the technique and proce¬ 


dure are presented inIPavlenko et al 


dwarf GJ 406. Eollowing iPavlenkol (1199 


in the case of the M6 
we computed par¬ 
tial pressures of some molecules and atoms which exceeded the 
pressures of the gas-dust phase transition and decreased their gas 
phase abundances to the corresponding equilibrium values. 

The bes t ht to the observed low-resolution spectrum of 
SDSS1256 (iBurgasser et al.ll2009l) was obtained for an effective 
temperature of 2600 K and a log(g) surface gravity of 5.0 (cgs). 
To get a better fit we reduced the absorption of CrH theoreti¬ 
cal bands by factor two. This might be caused by a deficit of 
Cr atoms in the atmosphere or some depletion of Cr into dust 
particles. Then, we increased the absorption of CaH bands by 
factor two, likely related to the well known enhancement of cal¬ 
cium an d other alpha-e lements in the atmospheres of metal-poor 
dwarfs (lMagainl|1987h . Moreover, we added some phenomeno¬ 
logical dust opacity i n the modelling of theor etical spectra, fol¬ 
lowing the scheme of IPavlenko et al.l (12007 ah . The wavelength- 
independent dust opacity of total an optical depth of t = 0.3 is 
located in a cloud-like layer at T/rmraii = 0.001, i.e above the 
photosphere of the dwarf. These phenomenologically imple¬ 
mented parameters provide a much improved ht to the spectrum 
of SDSS1256, reproducing the strong CaH and TiO absorption 


bands between 6300 and 7200A, as well as individual metal lines 
in this region (right-hand side panel in Eig. |2l. The most no¬ 
table discrepancy between the synthetic spectra and the observed 
spectrum below 6300A is likely caused by incompleteness of the 
CaH linelist at these wavelengths. The full spectral energy distri¬ 
bution of the SDSS1256 from the optical up to the mid-infrared 
is also well-reproduced by our synthetic spectra (Eig. 1^. These 
adjusting parameters should be considered in the framework of 
more accurate investigation of abundances and chemical equi¬ 
librium in the atmospheres of metal-dehcient L dwarfs because 
our parameters are actually not self-consistent with the tempera¬ 
ture/pressure prohle that generated the synthetic spectrum. 

To assess the reality of the Li I absorption feature in the spec¬ 
trum of SDSS1256, we examined the region around the Li I ab¬ 
sorption doublet by hxing the adjusting factor for CaH to the 
aforementioned values and setting [Li/H] = -1, -2, -3, and -4 
dex (corresponding to abundances of Li of 10, 100, 1000, and 
10000 less than solar, respectively). We note that the Li/H = -3.0 
and -4.0 models are nearly identical because Li opacity has 
fallen below some other opacity source (mainly CaH). Our ob- 


^ kurucz.harvard.edu 

* http://www.pages-perso-bertrand-plez.univ-montp2.fr/ 
^ http://vald.astro.univie.ac.at/ vald/php/vald.php 
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Fig. 5. Left: VLT/FORS2 obs erved spectrum of S PSS 1256 (black line). Overplotted in red is the synthetic sp ectrum from the 
NextGen atmosphere models of iBaraffe et alJ (Il998h generated with the DUST546 code of IPavlenko et al.l (12007 ah . which provides 
the best ht to the data (we show a low resolution version for plotting purposes). The synthetic spectra, convolved by the instrument 
prohle, have been generated for a metallicity of [M/H] = -2.0, gravity of logg = 5.0, Li abundance of logN(Li) = -1.0. Both spectra 
are normalised to the observed spectrum around 7050A. Major lies and features present in optical spectra of metal-poor dwarfs are 
shown as well. Right: Zoom on the region around the lithium doublet at 6707.8A. The observed spectrum is shown in black while 
synthetic spectra with different Li abundances are displayed in colour. 



Fig. 6. Full spectral energy distribution of SDSS1256 from the 
optical all the way to the mid-infrared. Overplotted in red is the 


synthe tic spectra from the atmosphere models ofjAllarde^^ 


( 20011) generated with the DUST546 code of iPavlenko et al 
(12007 al) . which produces a satisfactory ht. 


served absorption is weaker than the two most depleted abun¬ 
dances, suggesting that our tentative detection is unlikely asso¬ 
ciated with Li and instead comes from a band most likely created 
by CaH or TiO. We can set an upper limit on the Li abundance 
in SDSS1256 to [Li/H] = -3.0. We conclude that low Li abun¬ 
dances in subdwarfs will remain a challenge to detect (see also 
Section lTTt . 

6. Luminosities and effective temperatures 

We determined the luminosities and effective temperatures 
of both subdwarfs f ollowing the procedure described in 
iBurgasser et al.l (l2008l) . Previously publish ed low-resolution 
optical and near-inf r ared spectral data (iBurgasseil l2004t 
IBurgasser et al.l l2007l l2009h spanning 0.6-2.5 microns were 


combined and scaled to 2MASS and SDSS photometry. These 
were joined to segments of a T g/-f -2000K, logg - 5.0, solar- 
metallicity BT-Settl08 model (lAllard et al.l 120121) scaled to 
match SDSS and WISE photometry and Spitzer/IRAC measure¬ 
ments for 2MASS1626 0. The combined spectrum was then 
scaled to absolute Hu xes using the parallax measurements of 
ISchilbach et al.l (l2009l) . and integrated over 0.3-30 microns (ac¬ 
counting for 99.96% of the total spectral model Hux) to calcu¬ 
late bolometric huxes and luminosities. We obtained radius es ¬ 
timates from the ev olutio nary models of iBurrows et al.l (1200 ll) . 
iBaraffe et akl (l2003l) . and ISaumon & MarlevI (l2008h . assuming 
ages of 5-10 Gyr and the measured luminosities. Values are 
listed in Table |4] with uncertainties that include spectral, pho¬ 
tometric and astrometric errors, and scatter in model radii. We 
also experimented with different spectral model sets, model 
Te//s spanning 1700-2400 K, and subsolar metallicity models 
([M/H] = -l, -2); however, because the spectral model seg¬ 
ments encompassed a small portion of the spectrum these vari¬ 
ations were much smaller than the observational uncertainties. 
Indeed, varying the model parameters (i.e. the temperature by 
+500 K and the gravity by +0.5 dex) did not change the result¬ 
ing effective temperature by more than 15 K. 

Our analysis reveals that SDSS 1256 is 0.2 dex less lumi¬ 
nous and nearly 200 K cooler than 2MASS1626, despite be¬ 
ing classihed 0.5 subtype earlier. This is conhrmed by its 
fainter absolute infrared magnitudes (Mvri = 11.49+0.30 mag 
vs Mw 2 - 10.46+0.09 mag). The seeming reversal in temper¬ 
ature trend with spectral type may rehect differing metallici- 
ties between these ob jects or possibly unresolved multiplicity 
(iGoldman et al.ll2008l) 


® Measured values are [3.6] = 13.27+0.04 mag, [4.5] = 13.19+0.04 
mag, [5.8] = 13.03+0.06 mag, and [8.0] = 12.95+0.04 mag (DDT-225, 
PI Burgasser). 
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Table 4. Bolometric luminosities {\ogiQ{LtoilL q)), effective tern- 7x10® cm, and a disk passage time of 1 Myr {W - 200 pc/Myr). 
peratures (Tes), and radii (R) for SDSS1256 and 2MASS1626. For this rate: 


Name 


Teff 

R 


dex 

K 

Rg 

SDSS1256 (sdL3.5) 

-3.93±0.08 

2000±100 

0.090±0.003 

2MASS1626 (sdL4.0) 

-3.73±0.02 

2180±30 

0.095±0.002 


7. On the presence of lithium in halo brown dwarfs 

7.1. Weak Li I in L subdwarfs from accretion 


Given the absence of Li I absorption in excess of the overlapping 
molecular absorption feature, we infer that both L subdwarfs 
have depleted their lithium through core fusion an d have masses 
above 0.065 Mq. Indeed, evolutionary models of iBaraffe et al.l 
(Il997h show that metal-poor dwarfs with effective temperatures 
of 1900-2400 K have masses around 0.08 Mq (Tabled and have 
already depleted their lithium and lie on the main-sequence; the 
stellar/substellar boundary being around 0.06 Mq (with a spread 
of about 10%) at metallicities between -2.0 and -1.0 dex. 

However, it would be conceivable that a low abundance of 
Li could be maintained in the photosphere from either the inter- 
sellar medium (ISM) or circumstellar material. We can estimate 
the observable abundance in either case assuming a steady state 
photospheric composition, with accretion providing a Li source 
and convection and core fusion providing a Li sink: 


d ifLi.ph ^ph ) 
dt 


M, 


ph- 


dfhi.ph 
dt 


' fLi,a 


dAla, 


dt 


Mph 

■Mph-^^o 

^dep 


( 1 ) 


^Li.ph 
NLi,ISM 


10^ 


(4) 


For an Li ISM abundance consistent with primor dial Big Bang 
nucle o systhesis (log A (L/) = -t-2.6 or [Li/H] = 1.5; ISnergel et al.l 
120071: iLoddersI l2003h or hot halo sta re {Xo gNfLi) -+2.Q or 
[Li/H] =0.9: ISDite & SDit3ll982l : iFieldsIlTOllh this fraction cor¬ 
responds to photometric abundance of [Li/H] = -4 to -3, at or 
below our inferred confusion limit with CaH. Hence, it is un¬ 
likely that Li accreted from the ISM will be in sufficient abun¬ 
dance to be deteced in Li-burning L subdwarfs. 

In the case of circumstellar accretion. Solar meteoritic ma¬ 
terial has a Li abund ance 20 times higher than hot halo stars 
{\ogN{Li)isM - +2.2: lLoddersl200^ . so the corresponding mass 
accretion rate can be 20 times lower to achieve the same basal 
abundance. For a detection threshold of [Li/H] = -3, we can 
rule out circumstellar accretion onto SDSS1256 at a rate higher 
than dMcsmIdt » 6x10“^^ Mo/yr » 10® g/yr, roughly equiva¬ 
lent to the meteoric accretion rate of Earth. While not enough 
is known about the formation of solid bodies around these low 
mass, metal-poor stars to make a quantitative assessment, this 
limit suggests very little circumstellar material remains around 
these sources. Note that accretion of a larger planet could tem¬ 
porarily boost the lithium abund ance, a mechanism that has been 
proposed for Li-rich red gi ants (lAlexanderlll967Hlsraelian et al.l 
1200 It ICarlberg et al.ll2010h : however, such an event must have 
happened quite recently (<10 Myr) for an excess to be observed 
before the accreted Li is convected to the core and destroyed. 


Here, fu.ph is the steady state mass fraction of Li in the pho¬ 
tosphere; fu,acc the mass fraction of Li in the accreted mate¬ 
rial; Mph is the constant mass of the photosphere, estimated as 
10“*®M,o, ss 10 “"Mo (iBurrows & Lieberl]l993h : dMacddtis the 
instantaneous mass accretion rate in Mo/yr; and t^ep is the deple¬ 
tion time scale for Li by downward convection and core burning, 
assumed here to be 10 Myr. These assumptions yield 


fLi.ph NLi.ph tdep dM^cc nigtlMacc 

— - » - r ;- - — 10 — - — yr/AlQ 

jLi^acc N Li.acc f^ph dt dt 


( 2 ) 


The ISM accretion rate can be estimated from simple ballistic ac- 
cretionO i.e., a projectile of cross section nR^ plunging through 
the ISM at a relative speed v =» 300 km/fl In this case. 


dM, 


ISM 


dt 


rtHniHitR^v ss NHttiHrtR^Itp 


Mo/yr 


(3) 


where the first expression is the instantaneous accretion rate and 
the second is averaged ove r the disk passage, assu ming tin - 
1 cm^^ and Nh — 10^' cm“^ (iKalberla & KerTj2009t) . R^ R jup — 


’ Whil e ISM accret i on is commonly m odeled as Bondi-Hoyle accre¬ 
tion (e.g. nbshiiiri98ll : [Frebel et al.ll2009h . the dynamics of a low-mass 
halo subdwarf moving through the ISM at high speed does not satisfy 
the conditions for this kind of accretion, as the dynamic scale, GM/v^ 
S5 10'° cm, is much smaller than the ISM mean free path, ((Tjsmiiism) ^ 
i g 3x10*^ cm, assumi ng (Tism = n x 10“*® cm^ and hism = 1 cm“^ 
jBondi & Hovlein'9441) . Indeed, this condition is rarely satisfied for low 
mass stars. 

* The ISM is at rest relative to the local standard of r est, probably 
accurate to within -I-/-30 km/s (e.g. iFrisch & Slavinl[2006l) . so the rela¬ 
tive speed is just the absolute value of (U, V, W) = 325 km/s which we 
rounded down 


7.2. Searches for Lithium in Cooier Subdwarfs 


Given the lack of Li in mid-type L subdwarfs, can we con¬ 
firm the substellar nature o f halo brown dwarfs with the lithium 
test (iMa gazzuet ^ 199 li [Rebolo et al.l Il992t iMagazzu et al.l 


Il992r Basri et al. 1996!) in the near fut ure? For solar-metallicity 


stellar and substellar ch emistry (e.g. IBurrows &Sharpl 119991: 
iLodders & Feglevll2006l) . Li is predicted to be in atomic form 
above ~1500-1700 K at 1 bar pressure but is converted into 
molecules at lower temperatures. Hence, we would expect chem¬ 
ical depletion of Li for ultracool dwarfs below 1500 K, which 
seems supported by obs e rvations of the coolest L and T dwarfs 
jKiAgatricket^ [2008t iKing et aP 120101 : iFahertv et al.ll2014l 
iLodieu et al .11201 4). 

Theoretical models locate the hydrogen-burning limit at 
higher masses (~0.09 Mq for Z-0 vs 0.075 Mq for solar 
meta llicity) and warmer temperatures for sub-sola r metallici¬ 
ties (IBurrows et al.lll99^ IChabrier & Baraffjll997h . From our 
results for the two mid-L subdwarfs presented in this paper, 
we can conclude that the Li burning limit lies below 2000 K, 
equiva lent to masses below ~0.08 Mq. The low-metallicity mod¬ 
els of IBaraffe et al.l (Il997h predict a Li boundary minimum 
masses around 0.06 Mq and temperatures below 1700-1600K 
for subsolar metallicities and ages older than 1 Gyr, in the range 
for coolest L subdw arfs (e.g. 2MASS J05325346+8246465; 
iBurgasser et al.l l2008h . However, being substellar, the atmo¬ 
spheres of these objects become exceedingly cold at the older 
ages expected for halo dwarfs. An object with [M/H] = -1.0 
(-2.0) at the Li boundary mass has temperatures of <1000K 
(1050 K) for ages larger than 5 Gyr; i.e. it is a T subdwarf. If 
Li chemistry follows solar metallicity trends, the atomic form 
will likely be depleted in such objects. On the other hand, it 
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is possible that the Li condensation temperature is metallicity- 
dependent, just as mineral/metal c ondensate formation i s in¬ 
ferred to be from observables (e.g. iBurgasser et alJl2007h and 
models (e.g. IWitte et ^l2009l) . The search for Li absorption in 
even cooler subdwarfs is thus critical for constraining metallic- 
ity dependencies on both the Li-burning and hydrogen-burning 
mass limits, and atmospheric chemistry. 

8. Conclusions and outlook 

We presented optical spectra of two mid-L subdwarfs and set 
upper limits on the pseudo-equivalent widths of the Li I doublet 
at 6707.8A. We can summarise our results as follows: 

- we report the detection of a feature in the spectrum of 
SDSS1256 aligned with Li I at 6707.8A , but attributable to 
CaH. This sets a confusion limit for detectable Li abundance 
of [Li/H] = -3 

- we set an upper limit of 90 mA on the EW of the Li dou¬ 
blet for 2MASS1626, placing both objects in the Li-burning 
regime 

- we revise previous radial velocity measurements for both ob¬ 
jects and derive new space motions and Galactic orbits, the 
latter confirming both sources as members of the inner halo 
population 

- we find that ISM accretion of Li-rich material would provide 
a basal abundance below our confusion limits, and constrain 
circumstellar accretion to be less than 10® g/yr on SDSS1256 

- we reproduced the observed spectrum of SDSS 1256 based 
on NextGen atmosphere models with reduced metallicity 
and factors adjusting the amount of TiO and CaH molecules 

The next step is to attempt a detection in cooler L and T 
subdwarfs to reveal their true nature and place empirical con- 
straint(s) on low-metallicity evolutionary models. Further pho¬ 
tometric searches to uncover T subdwarfs will likely be needed 
in order to investigate the presence (or absence) of Li in metal- 
poor brown dwarfs. Moreover, additional work on the conden¬ 
sation of major chemistry features as a function of temperatures 
and pressur es in me t al-poo r atmospheres is needed, similar to 
the work of lLodde^ (1 19991) for solar-type brown dwarfs. 
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